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I. INTRODUCTION 


Contained in this report are the results of a study on the evalua- 
tion and reduction of noise impact to a community due to aircraft land- 
ing and takeoff operations. This work is a continuation of the methods 
and results of a previous study done by the same authors (under NASA 
Grant NSG-1509, reference 1). For completeness some repetition of the 
earlier work is included. 

The previous work considered only a single aircraft using a single 
approach/landing trajectory. Models of population distribution, air- 
craft noise signature, and aircraft flight path were developed, and a 
suitable annoyance model adopted. A performance index to be minimized 
was formed from the annoyance model and constraints. The current study 
has examined the case of multiple aircraft, flying on several trajector- 
ies, for either the case of approach/landings or for takeoffs. A su- 
perior, more realistic model of the flight path has also been developed. 
As in the earlier work, the annoyance criterion used is the noise impact 
index (Nil). The algorithm developed has been applied to Patrick Henry 
International Airport. 

Discussions of the various models, the performance index, optimiza- 
tion methods, and results appear in the following sections. 
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II. PROBLEM FORMULATION 

OVERVIEW 

The problem considered is that of determining the "best" set of 
aircraft landing and/or takeoff paths from any airport which minimizes 
the noise impact on the surrounding community. There are five major 
aspects of this problem which must be modelled: (1) aircraft noise 
signatures, (2) population distributions, (3) a cost function or per- 
formance index, (4) the aircraft flight paths, and (5) constraints on 
the aircraft (based upon aircraft dynamics), passenger comfort, safety, 
and maximum noise exposure for any population group. In addition, a 
flight path optimization scheme must be adopted. A modular concept has 
been employed so that any section of the problem may be modified with 
relative ease. The following sections describe each of these in detail. 

A. Aircraft Noise Signatures 

An aircraft noise signature gives a description of the noise ema- 
nating from an aircraft. Many such representations are available. The 
one adopted here is a simple model to facilitate computation; however, 
it can be replac d with more complex and accurate models. One such 
model is available through the use of the Aircraft Noise Source and 
Contour Estimation computer programs (see references 2, 3). The air- 
craft noise signature used in this study is obtained using data from 
reference 4. Here the effective perceived noise level (EPNdB) is given 
as a function of slant range to the closest point of approach for a 
variety of aircraft. A typical plot of the slant range variation for 
two different aircraft is shown in Figure 1. These data were fit using 
standard least squares techniques to yield an expression for EPNdB given 
by 
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Effective Perceived Noise Level (EPNdB) 




EPNdB = 115 - 22.5 x ( - lant ran8 ~- ---) 

500 

This equation is used for calculation of the maximum noise level at each 
location for a flyover. A typical footprint for a straight-in approach 
along a 3-degree glide sl^ce is shown in Figure 2. For other aircraft, 
similar experimental data must be sought in the literature. 

B. Population Distribution Model 

To model the distribution of population, a map of the community is 
overlaid with a grid and the population in each grid section is deter- 
mined. The population distribution within each section is assumed to be 
uniform. Several grid geometries were examined (see Figure 3). The 
geometries include: (l) rectangular sections of equal size, (2) rec- 

tangular sections whose dimensions increase with distance from the 
airport runway, and (3) concentric circles divided by several radial 
lines. The second scheme was chosen since it requires fewer rectangular 
sections than the first and is somewhat easier to implement than the 
third. Computer time required for determining the optimum trajectory 
varies directly with the number of grid sections. This results in the 
desire to minimize the number of blocks in the grid. Furthermore, since 
the noise levels decrease with distance from the aircraft and the air- 
craft has higher altitude when farther from the runway, the need for 
high resolution of the population density diminishes with distance from 
the airport. Grid blocks with larger area may then be used when farther 
away from the airport. 

Within a grid section, the population is determined by use of the 
SITE II system (reference 5), available on the CDC 7600 computer at the 
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Figure 2. Noise Footprint 
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NASA-Langley facility. This system requires as input the latitude and 
longitude of a reference point and the coordinates of the corners of 
each rectangular section. Although SITE II allows for simple retrieval 
of 1970 census data, there is some question about its resolution capa- 
bilities for small grid sections. In addition, in rapidly growing areas 
the population data may lag the actual population. The SITE II program 
is capable of producing detailed census information as shown in Figure 
4; however, for the present analysis only the population information is 
used, as indicated. 

C. Flight Path Mode l 

There are two ways in which the trajectory of the aircraft may be 
determined. In one, a discrete time integration of the equations of 
motion (with control deflections) yields point by point spatial coor- 
dinates and orientation. Although this allows the flexibility of ex- 
plicitly including control constraints as well as dynamic constraints 
(e.g. maximum roll angle), it requires that a considerable number of 
states of the system be stored in the optimization routine (i.e. each 
point of the trajectory in discrete form). In the multi-aircraft, 
multi-trajectory problem investigated here, such storage requirements 
are prohibitive. 

Thus, another method was adopted which uses only the functional 
form of the trajectory to describe the flight path. Two possibilities 
have been investigated: (1) a truncated Fourier series representation 
and (2) a scheme of line segments joined by smooth arcs. 

The Fourier series has the advantage of being able to represent any 
smooth function over a finite range reasonably well when the series is 
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truncated after a few terms. However, it is not able to represent 
functions with slope discontinuities without introducing "waviness" into 
the approximation. A large number of terms are needed to reduce this 
effect. The line segment representation does not have either of these 
features; however, It can approximate very well functions which describe 
the types of paths aircraft customarily fly. 

The first method begins by generating a starting path which goes 
from the initial trajectory point to the desired runway, ending up with 
the proper heading, i.e., the aircraft velocity vector is aligned with 
the runway. This starting trajectory is generated using the following 
equation: (see Figure 5) 

y s (x) = I« f (x-x p ) + (y p -y 0 )l <*P (- c, )1 + y 0 

For the vertical motion a simple three-degree descent path was assumed. 

Next, the first five Fourier sine harmonics are used to introduce 
deviations from the starting path. The coordinate system is scaled so 
that each of the sine functions contributes zero deviation at the end 
points. Therefore, if the starting path satisfies the boundary condi- 
tions, then the path with the deviations will also. An exponentially 
decaying factor is used to eliminate heading deviations at the final 
point. 

With the deviations, the equations for the path become 

y - {Z a. siu [in ( - ■ ■ - ) ] } [1 - exp ( + y (x) 

i=l f o 2 

N x ” x x * x 

z = {2 p. sin [in (——-)]} [1 - exp (— g — -)] + z„(x) 
i=l f o 2 
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where the ol and are the unknowns to be determined. 

The second flight path model represents the trajectory as a chain 
of line segments extending from the initial to the final point. Each 
corner between two segments is "smoothed" with a circular arc whose 
radius is large enough to insure that the aircraft can perform the turn 
(see Figure 6). The unknown variables to be determined are the coor- 
dinates of the line segment intersections (corner points). For the 
starting trajectory, the corner points lie equally spaced along a line 
through each pair of initial and final points. The number of line 
segments and hence, the number of corner points, per trajectory is 
determined before the optimization begins. This number is generally 
small (3 to 5) so that the pilot is not overburdened with required 
maneuvers . 

Both models of the flight path have the advantage of requiring only 
a small number of parameters to describe the trajectory. This reduces 
the optimization problem from a variational one to an ordinary one, but 
care must be taken to see that the various constraints in the problem 
are met. 

D. Constraints 

The use of a functional form of the flight path for the trajectory 
requires the reformulation of constraints into parameters which, can be 
used in the optimization. This is accomplished by translating the 
steady state solutions of the lateral and longitudinal perturbation 
equations into geometric constraints. For a detailed derivation of 
these, see the final report for 1979, Appendix A of reference 1. The 
constraints are incorporated by determining maximum curvature and slope 
parameters as a function of aerodynamic and physical constraints. 
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Similar expressions are given in the appendix (referred to above) for 
constraints on aileron, rudder, and elevator deflections, flight path 
angle and pitch rate limits. 

In addition to the aircraft constraints, there are passenger com- 
fort considerations (e.g. max bank angle), maximum noise exposure 
levels, and a minimum separation distance between multiple trajectories. 

All of the constraints are listed below: 

I. Aircraft Dynamic Constraint: 

d 2 y 
dx 2 


Lateral 


1 ♦ <£> 


C 1 + C 2 C 3 


v min (6r lf 6r„, 6rJ 

avg 


where V = average velocity of aircraft, Cj, C 2 , C 3 
are constants for a given aircraft, and the 6r's depend upon 
maximum bank angle and maximum rudder and aileron deflections 
for a given aircraft. 


Longitudinal: tan y 


max 


dz 

dx 


< tan Y d 


max 


where y r 
angles , 


and y . 


are the maximum climb and descent 


max 


max 


II. Passenger Comfort Constraint: 
3/2| 


a + <&>' i 

dx 


d 2 y 
dx 2 


V 2 

avg 

C 4* 


where C^ = 1.9 for 90% passenger satisfaction, 4.5 for 80% 
satisfaction, g = acceleration due to gravity, and V 
average velocity of aircraft during the turn. 


avg 
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III. Threshold Noise Constraint 

No populated area may receive noise in excess of 95 dB more 
than N percentage of times per day, where N is a fixed percentage 
of the number of flights per day. N is made as small as possible 
for any given case. 

IV. Minimum Separation Constraint 

A minimum distance of 800 meters (% mile) must be maintained 
between any two trajectories at all points (except very close to 
the runway, where all trajectories must converge). 

E. Cost Function 

A large number of criteria have been proposed to evaluate noise 
annoyance (e.g., EPNdB, NNJ, sleep interference index, speech interfer- 
ence index, etc.). The recent trend in noise assessment work is toward 
a universal measure — the noise impact index (Nil). This measure is a 
weighted day-night model which accounts for population density. It is 
described in detail in reference 6. Briefly, the total population 
exposed to each incremental average day-night model sound level is 
multiplied by the weighting function for that level. The weighting 
factor W(L^ n ) , multiplied by the population exposed to that L dn , is 
summed and normalized by the total population giving the Noise Impact 
Index for the area: 


Nil 


* P < L dn> W < L dn> 
L dn 


1 

L 


A plot of W(L dn ) 


P < L dn> 

dn 

appears in Figure 7. 
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The cost function or performance index for the optimization pro- 
cedure is taken to be the Nil plus penalties for violating constraints. 
Basically, the optimization procedure is set up to "drive" the aircraft 
trajectories to the path which will minimize the Nil and at the same 
time, not violate any constraints. As an example., the constraint of 
flight path angle not exceeding a maximum descent angle, Y d > nor a 
maximum climb angle, y c > is written as 

tan Y < 
c 

Each is converted to a penalty which is added to the Nil in the form 
Cost = Nil + KjPj + K 2 P 2 

P x = {max[0, (tan Y c “ dz/dx )]} 2 

P 2 = {max[0, ( dz/dx - tan Y d )]] 2 K^, K 2 = constants 

As is seen, for values of the flight path angle within the allowable 
range, no penalty is added; however, for values outside this range, the 
penalty and thus, the increase in cost, is great. Other penalty terms 
are added in a like manner. 


dz 

dx 


< tan Y d 
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III. OPTIMIZATION 


The optimum set of trajectories is determined by calculating values 
of the unknowns (the 0^ and in the Fourier series model, the corner 
points in the line segment model) which minimize the total cost (Nil 
plus penalties). Two optimization algorithms have been examined: the 
method of steepest descent and the Davidon-Fletcher-Powell method. An 
example of steepest descent is given below. Basically, the method 
computes the gradient of the cost function, C, with respect to the 
unknown parameters and then searches along the negative gradient direc- 
tion for values of the parameters which reduce the cost. 

In Figure 8, the point Lj represents the set of parameters which 
corresponds to the starting trajectory. The arrow points in the direc- 
tion of the negative gradient of C (i.e., the direction of decreasing 
Nil). Searching along this direction will yield a new point I^ which 
corresponds to a new trajectory with lower Nil. The process of comput- 
ing gradients and searching continues until the cost converges to within 
a specified tolerance. In this example, the sequence begins at and 
converges to L*, where the Nil is an absolute (or global) minimum. 

Consider, however, the case where the starting trajectory is char- 
acterized by the point Q^. The optimization process will converge to 
the point Q* , which is a relative (or local) minimum. The trajectory 
characterized by Q* does give a lower Nil than the starting path at Q^, 
but the Nil at Q* is still higher than that at L*. 

In this example, it is easily seen that if the starting point lies 
in Region I, convergence to the global minimum at L* is assured (like- 
wise for Region II and the convergence to the local minimum at Q*). The 
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only way to insure that the point I* is found is to execute the opti- 
mization algorithm a number of times with different starting points, as 
indicated by the open circles in Figure 8. There is the possibility 
that the cost function has a "sharp" global minimum, such as at the 
point R*. In such a case, it is likely that none of the starting points 
chosen would result in convergence to R*. From a practical point of 
view, though, it is not important that the true global minimum at R* is 
not found. The range of parameters defining the sharp "well" at R* is 
so narrow that a pilot could not deviate from the optimal path charac- 
terized by R* without greatly increasing the Nil. Simply stated, the 
only optimal path of interest is one whose resulting Nil is not overly 
sensitive to slight variations in the path. 

The steepest descent algorithm has the disadvantage of giving slow 
convergence near the optimal set of unknown parameters; however, signif- 
icant reduction in the cost (Nil) does occur during the first few itera- 
tions. A superior algorithm is the Davidon-Fletcher-Powell method, 
which gives good convergence near the optimum. This method has been 
employed in this study with satisfactory results. A detailed descrip- 
tion of both optimization methods appears in reference 7. 

A. The Optimization Algorithm 

A computer code has been developed which implements either of the 
optimization methods described above. Figure 9 shows a flow chart for 
this code. Initial data (population map, aircraft constraints, initial 
and final aircraft positions, etc.) are required for each configuration 
of trajectories and aircraft at a given airport. An initial set of 
trajectories is either supplied by the user, or a default set is gen- 
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Figure 9. Flow Chart of the Flight Path Optimization Algorithm 
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erated by the program. The optimization then begins, with successive 
values of the cost being compared after each iteration. When the dif- 
ference between successive values is less than a defined stopping cri- 
terion, the process terminates. 

The code has been written in modular form so that any of the vari- 
ous models (population distribution, cost function, etc.) may be up- 
graded or modified easily without making major changes in the code. As 
an example, the noise impact in each population section requires the 
computation of an integral. While this integral is usually approxi- 
mated, a more accurate calculation can be made with the simple addition 
of a subroutine to the program. 

Appendix A '":ontiiins the FORTRAN code as written for a CDC Cyber 172 
machine. 

B. Results 

All of the cases discussed here involve the Patrick Henry Interna- 
tional Airport in Hampton, Virginia. The SITE II program was used to 
generate the population data for each block as shown in Figure 10. The 
three entry points referred to, Swing, Franklin, and Cape Charles, are 
the check points indicated on the ILS approach plate (figure 11). 

Reduction in the Nil at Patrick Henry Airport is limited by the 
population distribution. As indicated in Figure 12, most of the people 
are located in blocks near the runway. During takeoffs and landings, 
these people will be affected by aircraft noise regardless of the tra- 
jectories flown. 

(1) The Swing and Franklin entry points are used simultaneously for 
approach/landings. With the Fourier series model of the flight path, 
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Sine Harmonic Representation of Flight Paths 




the results obtained are shown in Figure 12. Details appear in Table I. 
As is easily seen, there is an unnecessary amount of waviness in the 
trajectories far from the runway. This is caused by the fact that the 
Fourier series is truncated after five terms. More terms could be 
included but more computation time would be required. Thus, the line 
segment model of the flight path has been adopted and is used in all the 
following cases. 

(2) A single trajectory, with one Boeing 707 flying, is determined 
using the line segment model. The results are shown in Figures 13a and 
b. Both the Swing and Franklin stations have been used as entry points. 
There are three segments in each trajectory, requiring only three turn- 
ing maneuvers from the pilot. This is clearly more realistic than the 
type of path produced by the Fourier series model. A comparison of the 
results of the two models shows that the line segment scheme yields 
slightly higher Nil values (3-5% higher than in the Fourier series 
representation); however, the Nil is reduced, compared to existing 
approach paths, by 4-6%. 

(3) Multiple aircraft on multiple trajectories are investigated. 
Figures 14, 15 and 16 show the results for two, three, and four segments 
per trajectory. The reduction in Nil ranges from 4 to 5%. Details 
appear in Table II. 

(4) The multiple aircraft, multiple trajectory case is repeated 
(with three segments per trajectory) using Gaussian quadrature to eval- 
uate the integral in the Nil computation; a 6% reduction is seen. 
Figure 17 shows the difference that results when this integration is 
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| Figure 15. Three Segments per Trajectory (Landing) 
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Figure 16. Four Segments per Trajectory (Landing) 
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Table IX Line Segment Representation of Flight Path 


(II) Lina Begment representation Of flight path 


(A) Balanced dlatrlbutlon of aircraft on each trajectory* 
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of 
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trajectory 

Annoyance 
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REMARK 


Swing entry 
traejetory 

Franklin entry 
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entry trajectory 


igufe 

umber 

i 

3 
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(800, 0) 



Single A/C, 
single trajectory 

13a 
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Single A/C, 
single trajectory 
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2 
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threshold noise & 
trajectories sepa- 
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14a 
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14b 

3 
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Optimal trajectories 
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traj) using centroid 
approximation in 
Nil calculation. 

15 

3 
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Figure 17. Comparison of Trajectories Using Gaussian Quadrature 
(dotted line) and Centroid Approximation (solid line) 
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performed as compared to the customary approximation.* Details appear 
in Table II. 

(5) The aircraft mix on the Swing and Franklin trajectories is 
altered to be unbalanced (details in Table III). Figures 18a and b show 
that between the two cases of unbalanced aircraft distribution, the 
resulting change in the trajectories is slight and the change in NIX is 
only 0.07%. This may point to the existence of optimal "corridors" 
which are independent of the aircraft distribution. 

(6) All three entry points are used simultaneously for multiple 
aircraft, as shown in Figure 19. The important result here is that the 
optimum trajectory from Cape Charles is found to pass over !;he water, as 
should be expected. Details appear in Table II. 

(7) Some preliminary work has been done on the takeoff problem, 
For each takeoff trajectory, the end of the runway becomes the initial 
point, and the final point (approximately 30 km away) may be placed 
anywhere. Two final points in the region northwest of the airport and 
two in the southwest region were chosen. Optimal paths were computed 
for the different pair combinations (three segments per trajectory). 
These are shown in Figures 20, 21, and 22 with details given in Table 
IV. The pair giving the lowest Nil is shown in Figure 21. (Nil = 
1.425). 

(8) To help guarantee that the optimal set of trajectories is found 
by the searching algorithm, a method called "selective search" has been 
devised. Figure 23 shows a simple version of it. Basically, a number 

-•■Referred to as the centroid approximation, since the in a given 

population block is calculated at the centroid of the block and 
assumed constant over the entire block. 
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of trajectories are evaluated as candidates for the starting set input 
to the algorithm. The process consists of choosing several segments in 
different regions of the near-terminal area and evaluating the Nil for 
various combinations of them. The set of trajectories with the lowest 
Nil is taken as the starting set for the optimization routine. This 
way, it is more likely that the optimal set will be found. Such a 
selective search corresponds loosely with the choosing of different 
starting points in the example given earlier (figure 8) . 

CONCLUSION 

A method has been formulated which optimizes aircraft paths during 
approach or takeoff. Multiple aircraft flying on several trajectories 
can be considered. Models have been developed using available data 
where possible for the population distribution, aircraft noise signa- 
tures, noise impact, constraints, and flight path. An algorithm which 
uses either the steepest descent method or the Davidon-Fletcher-Powell 
method for optimization has been implemented and tested. The algorithm 
can 

1) Evaluate the noise impact of existing flight paths, 

2) Evaluate the noise impact of proposed flight paths, and 

3) Optimize the flight paths to minimize the noise impact, sub- 
ject to constraints. 

The method has been applied to the Patrick Henry International 
Airport area. Existing paths have been evaluated for noise impact and 
optimal paths were determined using either two or three of the terminal 
area entry points. Approximately 4.5% improvement in Nil was achieved 
over that of the paths presently used. The population is concentrated 
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near the end of the runway, though, and it is felt that even more im- 
provement in the Nil could be achieved at other airports. 
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PNALTY • PNALTY + PWEITl*((APtAXl(0..( AL0NG2-A) I ) **2*( APAXKO.. I A 1220 
1 L0NG2-B ) ) )**2) 1210 

IP (J.NE.2) PNALTY • PNALTY«PVEIT1»( AHAXl (0.. (ALQNG1* ALDNG2-AI 1240 
1 1 • • 2 1 1230 

C 1260 

C 1270 

C . . 1260 

C . CALCULATE TANGENTIAL POINTS. CENTER OF RADIUS ft AT J-TH CORNER . 1290 

C . • llt'O 

C 1119 

C 1120 

XTld.JI • X H < I. J I* AL0NG2* (XH(l*J“l)”XNd»J>)/A 1110 

YUd.JI • YHII* J |4AL0NG2»( YH( l» J-l )-Yhd> J ) ) /A 1140 

XT2(l#JI • XMd,J).AlONG2*(XNd»J*l)-XN(I»J> W6 1130 

YT2(I.J> • YNd.J>mCNG2*IYN(l#J»i>-YMd»J))/6 1160 

SfcG > S0RT((XT2d.J-l>-XTld.JI)**2«(YT2d.J-l>-YTld.J>>**2) 1170 

DISN2 • AL0NG2*C0S ( OCT A/2. 1 1180 

01SC2 • RAD I US/ SIN (BET A/2. I 1190 

X > (YTHI. J>*K12(I, JII/2. 1400 

Y • (YTl<l.J)tYT2(I.JI)/2. 1410 

XCENTKI I. J-l) • XM(I,J)*DI5C2*(X-XHa» JM/0ISN2 1420 

Yi.ENTRII.J-1) • YM(I.J).DISC2*(Y-YNd»J))/0ISH2 1410 

AL0NG1 • AL0NG2 1440 

0 • SQRT((XTld»J)-XT2d»JI)**2*(YTld»J)-YT2d»J))**2l 3430 

ANGLE ( I .J-l ) • AACOSIR AGIOS. RADI US. 0) 1460 

SEG • S£G*ANGL£(I.J-1)*RADIUS 1470 

3u OlS(l) • DISdltSEG 1480 

40 CONTINUE 3490 

A • S0RT(<XN<I»NScG)-XMd,NSEGU)>**2«(YMd.NSEGI-YM(I»NSEGUM»3500 
l *2) 3310 

B • SORTI (XH(l.NSEGAl)-XMI.NSEG42l)**2«(YN(I.NSEG*ll-YNII.NSEGA3320 
1 2) 1 • ♦ 2 > 1530 

C • SORT) (XPd,NSEGI-XNC»NSEG*2) I »»2A< YNII. NSEG l-YN (I»NSEG*2) I *1540 
1 *2) 3550 

ANG • AACOS(A.B.C) 3960 

PNALTY • PNALTYAPR£IT1*(ANAX1(G.»(40D.*AL0NG1-A) l»*2 ) *PWEI T2*( AH3570 
1 AXllt..(2.3e0579899-ANG)l**2l 3560 

DIS(I) • DISd)*S0RT((XT2d.NSEG)-XNd.NSEG*li I **2A ( YT2 ( I.NSEG 1-1390 
1 YNd.NSEG*l))**2l 1600 

50 CONTINUE 1610 

C 3620 

C 3630 

C . • 3640 

C . RETURN IF DYNANIC CONSTRAINT NOT SATISFIED . 1650 

C . • 1660 

C 3670 

C 3680 

IF (PNALTY. E3. 0. I GO TO 60 1690 
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I 

I 

I 


SU6R0UT I He COST 


71/172 IS 


FTN 4.7*415 


60/04/24. 10.20.16 


TOTAL • PNALTY 
RETURN 


TARE SAMPLES ALONG THE TRAJECTORY 


60 00 120 1 • l.NTRAJ 

OELZPM « (20(1)-2FI1I)/D1S(1) 

K • 1 

posit ( i.i.l) ■ xnii.i) 

POSIT ( It If 2 ) • YH(Ifl) 

POSIT ( If if 1) • 20(1) 

00 luO J • 2.NSEG 

DELX • (XT2 ( If J-l)-XTl(IfJ))/NSAMP 
DELY • ( VT2 ( If J-ll-YTll!f JII/NSAHP 
DO 70 L ■ XfNSAHP 

POST T ( If K*lf 1 1 • POSlTUfKfll-OELX 
POSIT ( If K* If 2 1 • P0SIT(IfKf2l-CELY 

POSIT ( If K«lf 1 ) • PUSimfK>ll-0ELZPn*S0RT(0ELX**2*DELY**2 
K • Ml 

70 CONTINUE 

IF (ANGLEdf J-D.E9.PI) GO TO ICO 


SAMPLE THE CORRECT ARC ON THE J-l TH CIRCLE 


ALFA1 • ATAN2UYT1I If J)-YCENTR(If J'lllt IXTKIf J)-XCENTR(W J 
)) 

(PLUS • 1 

AA • ALFA1* ANGLE ( If J-l ) /4. 

Xi 2 • XCENTR : U J- 1 1 *RAOIUS*COS ( AAI 
Y22 • YCENTR'5. J-l )*RAD1U5*S1N(AAI 
AOIS • S0RTUX22-XfUIf J I ) «*2* ( Y22-VHI If J)i**2) 

601S • SORTUXTKlf J >-XM( I t J I > **2*< YT1 ( If J 1-YHC I f J))**2> 

IF IBC1S.GT.ADIS) (PLUS • -1 
II • 0 

A • ANGLE ( It J-1 ) 

A • A-. 1745129252 
IF (A.LE.O. ) GO TO 90 
II • 11*1 

AA • ALFA1*1PLUS*1I*. 1/45329252 
FOSITUfMlt 1) • XCENTRIlf J-1)*RA0IUS*C0S(-AA> 

POSIT ( If K*lt 2 ) • YCkNTRdf J-1)-RA0IUS*SIN(-AA) 

0 • RADIUS*. 1745129252 

PCSlT(IfK*lf 31 • POSITU.Kf 3I-DELZPM«D 

K • Ml 

GO TO 80 

POSlT(lfK*l» 1J • XT2 ( I . J ) 

POSITUtRflf 2) • YT2 ( I f J ) 

D • RAOIUS»UNGLE(If J-l)-FLQAT(II)»fc. 1745129252) 
PCSlT(lfK*lfl) ■ P0S1T(I.K.1)-DELZPP«D 


1700 
1710 
1720 
.. 3710 
. 1740 
. 1750 
. 1760 
.. 1770 
1760 
1790 
1600 
1610 
1820 
1610 
1640 
1690 
1860 
1670 
1860 
1690 
1900 
) 1910 

1920 
1910 
1940 
1950 
.. 1960 
. 1970 
• I960 
. 1990 
.. 4000 
4010 
-1)4020 
4010 
4040 
4050 
4060 
4070 
4060 
4090 
4100 
4110 
4120 
4110 
4140 
4150 
4160 
4170 
4180 
. 4190 
4200 
4210 
4220 
4210 
4240 
4250 
4260 
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PAGE K 

OF POOR QUALJTY 


SUBROUTINE COST 73/172 TS 


UN 4.74485 10/09/24. 10.20.19 


\ 


I 

I 

I 

I 

I 

I 

I 





K • K*1 

4270 



US 

CONTINUE 

4210 




DEU • (XT21I»NSEGI-XN(I.NSEG*1)>/N,\AHP 

4290 

179 



CELT • ( V T 2 1 I.NSLG l-YHI 1 »NSEG*L 1 l/NSAHP 

4300 




00 110 L ■ l.NSAHP 

4310 




POSIT ( 1.K41.1) • POSIT tl*K*l»-0ELX 

4120 




POSIT ( I . K41.2 1 • POSITII.K.2I-0ELY 

4130 




POSIT ( 1. Ml. 3 1 • POSIT! I.K. 1I-DELZPPP$0RT (DEL X**2*DELY*92 1 

4340 

180 



K • K*l 

4190 



110 

CONTINUE 

4160 




NPOSlTdl ■ K 

4170 



120 

CONTINUE 

41(0 




IF (NTRAJ.EQ.il GO TO 270 

4190 

199 



DO 13U 1 ■ l.NHAP 

4400 




00 130 J • 4.5 

4410 




DO 130 K > l.E 

4420 




LQCaLU.KI • 0. 

4410 




ARRAYU.J) > LOCAL ( 1. Kl 

4440 

190 


130 

CONTINUE 

4490 




DO ill 1 • l.NTRAJ 

4460 




00 140 K • l.NHAP 

4470 




LOCAL (K. 2 1 • 0. 

4460 




LOCAL(K.l) • LQCAKK.2I 

4490 

199 


I,4u 

CONTINUE 

4500 


C 



4510 


c 



00 4920 


C 

• 


. 4930 


c 

• 

CALCULATE ANNOYANCE LEVEL 

. 4540 

200 

c 

0 

ANY BLOCK NITH ZERO POPULATION IS BYPASSED 

. 4950 


c 

• 

IN /LOCAL*' C0L1 • HIGHEST 707 NOISE 

. 4960 


c 

• 

C0L2 • HIGHEST 727 NOISE 

. 4570 


c 

• 


. 4580 


c 



to 4990 

209 

c 



4600 




NPOSITI • NPQSITIII 

4610 




DO 160 J • 1. NPOSITI 

4620 




DO 160 K • l.NHAP 

4630 




IP (ARRAY(K.3I.EQ.0.I GO TO 160 

4640 

210 



RANGE • SORT (( PCS IT 1 1. J. 1 l-ARRAY! K. 11 1*»2 ♦(POSIT ( I. J. 2 (-ARRAY 14650 




L K»2)I»«2AP0SIT(I.J. 31*421 

4660 




AL707 • 129.-29.*ALQG10(3.261*RANGE/200.) 

4670 




AL727 • U4.-22.5»ALOG10(3.2814RANGE/500. 1 

4680 




IP (AL707.LE.L0CAL(K.1)> GO TO 190 

4690 

219 



LOCALlK.ll • 11707 

4700 



190 

IP (AL727.LE.L0CALIK.2II GO TO 160 

4710 




LOCAL ( K. 2 I • AL727 

4720 



190 

CONTINUE 

4730 




00 19U K • l.NHAP 

4740 

220 

c 



4750 





.4 4760 


c 

• 


. 4770 


c 

0 

IN .LOCAL*' COL 3-9 • HO, OF OCCURANCES HIGHER THAN ALHAX 

. 4780 


c 

0 

DUE TO IRAJ. NUMBERS 1-3 RESPECTIVELY 

. 4790 

229 

c 

0 

COL 6-8 • TOTAL VIOLATING NOISE - ALHAX FOR TRAJ. 

. 4600 


c 

• 

1-3 RcSPECTIVELY 

. 4810 


c 



. 4 820 


c 



... 4830 


I 


l 

I 

I 
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1 

i 




subroutine cost 73/172 is 


FIN 4. 7,4(3 (0/04/24. 10.20. IS 


■ 


r 


i 

,4 




I 

I 

I 

I 


1 

I 


211 


211 


240 


24S 


2So 


255 


2VJ 


2o5 


2 70 


275 


140 


2 35 


C 4(40 

IF (ARRAYIK, 31. EO.J. I 00 TO 190 4(30 

IF (4ftfUm.6l.hE.il 00 TO 190 4(60 

IF ( 10C6K K. ll.Lc.ALHAX) 00 TO 173 4(70 

10C4L(K»I»2) • NP7C71 1 1 ,NN707I 1 1 4*60 

LGCAL(K,I«5I ■ ( LOCAL ( K. 1 )»ALNAX I *( N0767I I )*NN707( 111 4(90 

170 IF (LOC ALIK. 21. IE. ALMA XI 00 TO 1(0 4900 

LOCAL (K» 1.2 I • LOCAL (K»I,2l,ND727!ll ♦NN72 7 ( 1 1 4910 

LQC ALIK. 1,51 • L0CAL1K»I»5),(LQCAHK,2I-ALHAXI*IND727(II,NN72T4920 
1 (111 4930 

IOC IF (LQCAL(K,I,2).LE • FI ITHAX ) 00 TO 190 4940 

THkESH • THRE SH,TNE1T1*( AN AX1I0. .(LOCAL! K, 1,2 )~FL1TNAXI l**2 I ,T49S0 
1 WE1I2*IARAX1I0.,L0CAI(K,I,5||»»2I 4960 

190 CONTINUE 4970 

2C0 CONTINUE 49*0 

00 2 it, I > l.NTRAJ 4990 

PEOPLE • 0. 3000 

00 21U K • l.NNAP 3010 

IF (ARRAYIK, 31. EO.O. I GO TO 210 3020 

ARRAYIK, 4 I • ARRAY IK, 4), NO 70 7111*10.*, I LOCAL (K.l 1/10. I*NN707( I 1*3030 

1 10.** (1., LQC ALIK, 1 I /10. 1 AND727I 1 1 * 10. ** I LCC AM K, 2 1 / 10. I*NN727 ( 1 1 3040 

2 ♦ iU.»*(l.+L£ICAl(K»2l/10.l 3030 

PEOPLE • PE0PLE,ARRAV(K,3I 3060 

210 CONTINUE 3070 

C 30*0 

C 3090 

C . . 3100 

C . CHECK THRESHOLD NOISE CONSTRAINTS . 3110 

C . . 3120 

C 3130 

C 3140 

00 22. K • l.NNAP 3130 

IF IARRAY(K.3I.E0.3.I GO TO 220 3160 

AVN • 1U.*AL0G10<ARRAY1K>4)/NPLANE) 5170 

ARRAYIK. 4 I • AVN 31*0 

IF (AVN. LT. 20. .OR. ARRAYIK, 31. EB.O. I GO TO 220 5190 

ARRAYIK, 5 I • 3 . 36e-6»10.»*( .103 «AVN I /( .2*10. ** l.03*AVN I *1.431-4*3200 
1 1>.»*(.U6*AVNI I 5210 

ARRAYIK, SI • ARRAYIK, 3I*ARRAY(K»5I /PEOPLE 3220 

AN1 1 • ANII, ARRAYIK, 51 3230 

220 CONTINUE 3240 

C 3230 

C 3260 

C . . 3270 

C . CricCK FOR THRESHOLD NOISE VIOLATION AT EACH BLOCK DUE . 32*0 

C . TO CURB INAT ION EFFECT OF HULT1PLE FLYOVER • 3290 

C . .3300 

C 3310 

C 3320 

00 26C 1 • l.NNAP 3330 

IF (ARRAY I I.6I.NE.1. I 00 TO 26w 3340 

VIOLA • Cl. 5330 

DO 130 J • l.NTRAJ 3360 

VIOLA • VIOLA, LCCALIl,J»2> 3370 

230 CONTINUE 53*0 

IF (VIOLA. Le.FLITNAXI GO TO 260 3390 

00 24vl J • l.NTRAJ 3400 


I 

I 


I 


OFP^i 1 PAGB IS 

OF Po °R quality 


SUBROUTINE COST 


73/172 TS 


FTN 4. 7*4*5 


•0/04/24. 10.29.13 


240 


295 


300 


305 


31J 


315 


32* 


325 


330 


335 


34u 


IF (LGCAl(I.J*2).GT .FL1THAX I SO TO 260 
240 CONTINUE 

SUN • 0. 

00 23u J • l.NTRAJ 

SUN • SUHfLOCALI 1.3*3) 

250 CONTINUE 

THRESH • THRESH* TOE I Tl*( AHAX1 (0.. ( VIOLA-F LITHAX 
1 AX1(0..SUH)**2I 

260 CONTINUE 


CHECK SEPARATING CONSTRAINTS 


270 

ie j 


3410 

3420 

3430 

3440 

3430 

3460 

) ) **2 ) *TN£IT2*( AH5470 
54 BO 
5490 
3500 
3310 
3320 
3330 
3540 
5330 
5360 
3370 
3360 
3590 
3600 
3610 
3620 
5630 
3640 
3630 


290 


300 


IF (NTRAJ.NE.il CALI CRDSBVR INFOS IT. POSIT. CLOSE > 

SO TO 32C 
DU 26i 1 • l.NHAF 
i'O 2dC J • 1.2 
13CALU.J) • 0. 

ClNTINUe 
OU 2VU I • 1 . K 
DO 290 J • i.NMAP 
If (AkRAYIJ.JI.EQ.U.I GO TO 290 
RANGE > S0RT((P0SIT(1.1.1)~ARRAY(J.1))**2*(PQSIT (1.1.21-ARRAY ( J. 3660 
i ,!>)**2*PUSIT(1.I.3)»*2I 3670 

AL707 • 129.-25. *ALCG1U(3.281*RANGE/2C0. I 3680 

AL/2.' • 115.-22. 5* ALOGlU < 3. 2 61* RANGE / 5C0. ) 3690 

LOCAL ( J > 1 ) • AnAKKAL707.L0CALU.lt) 3700 

LOCAL ( J. 2 1 • A PAX 1 ( AL72 7. LOCAL ( J.2 ) ) 3710 

CONTINUE 5720 

PEOPLE • 0. 3730 

00 300 I * l.NHAP 5740 

IF (ARRAYII.3I.E0.0.) GO TO 303 3730 

ARkA V ( 1. 4 ) • ARRAY (I.4t«N07L7(ll*10.**(LOCAL(I.l)/10. )*NN707(1I*5760 

1 1C. **(l.*LOCAL(l.l)/lo.)*ND727(l)*10.**(LOC ALII. 2)/lC.I*NN 727(1)5770 

! *10. **(1.*L0C ALII. 21/10.1 5780 

PeOPLc • PE0PLE*ARRAV(I.3) 3790 

CONI INUE 3800 

UQ 310 I • l.NHAP 5810 

IF (AKKAY(I>3).£S.u.) GO TO 310 3820 

AVN • 1*.*AL0G10(ARRAY( I.4I/NPLANE) 3830 

ARRAY! 1.4) • AVN 5840 

IF (AVN.LT.55.) GO TO 31k. 3850 

ARRAY! 1.5) • 3. 3 6E-6* 10. ♦♦ ( . 103*AVN ) / ( . 2*1C . *• ( . 03 *AVN)*1. 436-4*5860 


310 

320 


iO.**(.Ce*AVN)l 

5870 

ARKaYII.5) • ARRAYII.3)*AKRAY(I.5)/PE0PLE 

5880 

AN1I • ANI I*ARRAY (1.5) 

3890 

CONTINUE 

5900 

TUTAL > ANII»PNALTY*THRESH*CLCSE 

5910 

IF (1GRAD.LC.1) RETLRN 

5920 

DO 350 1 • l.NTRAJ 

5930 

SXTK1.1) • XTK 1.11 

5940 

SXUII.il • XT2li.ll 

5950 

iYHIl.ll • YTKI.l) 

5960 

SYT2I1.1) • YT2II.il 

5970 


57 


SU6R0UT1NE COST 71/172 T5 


FTN 4.7*419 10/04/24. 10.20.19 


DO 134 J • 2.NSEG 9010 

SXCcNTRU.l-ll • XCENTR U. <1-1 I 9990 

149 SYCtNTRU.1-1) * TCENTRU.l-l) *000 

SXT1U.1I ■ KUIliJI *010 

SYUU.ll • YTIU.1I 6020 

SXT2U.il ■ XT2U.1I 6010 

SY12U.lt * YT2U.1I 6040 

190 S4NGl.cU.l-ll • ANGLE ( 1 .1-1 1 6090 

130 CONTINUE 6060 

NPOSITSUI • NPOSXTU) 6070 

IPUST • NP0S1TSU I 6060 

00 14w 1 • 1.IP0ST 6090 

199 SPOSIIU.l.ll • POSlTUil.il 6100 

SPUSITU.l. 2 I • POS ITU. 1.21 6110 

SPOSIIU.l.ll • POSIT U.l. 31 6120 

340 CONTINUE 6130 

190 CONTINUE 6140 

36il RETURN 6190 

END 6160 


49UJ0B CM STORAGE USED 


10.132 SECONDS 


SUBROUTINE C3ST1 73/172 IS FTN A.7M69 B0/06/2A. 10.20.19 


SUBROUTINE COST1 ( IGRAO#N#F.X#ANII»PNALTT#ClOSE»THRESHI 6170 

COMMON NTRAJ#NHAP#NiEG#XM(3»6l.7MI3#6)#ARRAVt576#9l#SP0SlTI3#i00#36le0 
1)»1C(3I#2F(3I»NP0SITSI1I 6190 

DIMENSION 1 IN# 1 1 6200 

5 00 10 1 • 1#NTRAJ 6210 

NSEG1 ■ NSEG-1 6220 

00 lo J • 1»NSEG1 6230 

00 10 K > 1#2 6269 

1 • U-l)»2«lNSEG-ll*(J-l)*2»K 6290 

10 IF IK.EO.ll XMIl#J»l) • m#l) 6260 

IP IK.E0.2I TH(I#J*1I • XIL.1) 6270 

1C CONTINUE 6260 

CPU COST (IGRAD»F»ANII*PNAt.TY.CL3SE»THRESHI 6290 

RETURN 6300 

19 ENO 6310 


61UOOB CM STORAGE USEO .167 SECONOS 


59 


F 





FUNCTION AACOS 71/172 TS 


FTN 6. 7*6(5 (0/06/26. 10.29.1S 


I 

I 

I 


FUNCTION AACOS U.I.CI 6120 

X • IA»»2*B»*2-C*»2>/<2.»A*B» 6110 

IF ( ABB 1 1 ). IT. 1. 1 1 GO TO 1C i>160 

WRITc <6*9010) X.A.B.C 61S0 

S STOP 6160 

10 IF (X. ST. 1.1 X • 1. 61)0 

IF (X.LT.-l.) X • -1. 61(0 

AACOS • ACOStXl 6190 

RETURN 6600 

10 C 6610 

9910 FORMAT I29H TROliBLX IN AACOS* X.A.B.C * *61 1PE16. 9.1X1 1 6620 

END 6610 


6101o( CM STORAGE USED .995 SECONDS 



] 

I 

I 

I 

I 

I 


I 

I 

I 

I 


| 

I 




60 


SUIR3UIINE HONIT 13/172 TS 


flN 4.7*4 a5 10/04/24. 10.20.19 


lb 


19 


20 


29 


30 


39 


SUBROUTINE HONIT < I T. X.N. F. ANI I .PNALTY. CLOSE. THRESH) 4440 
CQHHGN NTR 4 J.HH 4 P.NSEG.XH 1 3.6 1 >TH( 1.61. ARRAY (9T6.9l.Sr0Sim.100. 36490 
l).ZO(3).2F(3 ) .NFOS ITS ( 3 1 6460 
LOHHON /PRINT/ SXCENTRC1*3)*SYCENTR(3»3)#SXTI(3.4).SYT1( 1. 4 ) » SXT2(6470 


13.4).SYT2(3.4).SANGLE(3.3) 

01HENS1QN X(N.l) 

00 lb 1 " l.NTRAJ 
HSEG1 • NSEC-1 
00 1C. J > 1.NSEG1 
00 10 K • 1.2 

L > ll-l)*2*(NSEG-ll*(J-ll«2*K 
IF (K.EO.l) XH(I.J.l) • XU.i; 

IF (K.E0.2I YH( I. J*1 I ■ XCl.lt 
10 CONTINUE 

WRITE (6.9G19) IT.F.4HU.PHALTY. CLOSE. THRESH 
DO 3u I • l.NTRAJ 

WRITE (6.9020) I »XH( 1 . 1 ) . YHC 1.1 1 
PII ■ i./ATAN(l.)/4. 

DO 2b J • 2.NSEG 


4410 

6490 

4900 

6910 

6920 

6910 

6940 

6990 

6960 

6970 

6960 

699C 

6600 

6610 

6620 


20 


30 


WRITE (6.9030) XHU. J I »YH( I. J ) .SXCENTR< 1. J-l t.SYCENTRC I. J-l ). S6630 
XT1(I»J)»SYTKI» J).SXT2(I.J).SYT2(I.J)»SANGLE(1. J-1I*PII*160. 6640 
CONTINUE *690 

WRITE (6.904b) XH( 1 .NSEG.l) > YH( I.NSEG«1 ) 6660 

NP0S1TI ■ NPQS'TSCII *670 

WRITE 16.90901 ( ( SPOS IT ( I. J.K ) >K>1. 3) « J«1.NP0S1T I ) 6660 

CONTINUE *690 

RETURN 6700 

* 6 7 1C 

9013 FORPAT (1H1.9X.11HITERATI0NI .12. /.12X.12H TOTAL COSTI . 1PE16.9. /. 16720 
12X. 17HANN0YANCE (NI1H , 1PE16.9. / > 12X. 31HPENALTY ON DYNAHIC C0NSTR6T30 
2A1NTI . 1Pc16.9./.12X.34HPENALTY ON SEPARATING CONSTRAINT! .1PE16. 96740 
3. /. L2X.26HPEMLTY ON THRESHOLD NOISE! .1PE16.9./I 6790 

9J2J F0RHA1 (12X.16HFLIGHT PATH NOi . I 1. /> 14X. 10HC0RNER PT.» 14X. 16HCENT6760 
ItR OF CIRCLE. 6X.19HTANGENT1AL PTS. .33X. 13HANGLE (DEGREE I • /.14X. 1H( >6770 
2FB.1.1H..FE.1.1HI) 6760 

9030 FORHAT ( 14X, 4 ( 1H (» F8. 1, 1H. ,F6.1. 1H). 9X1 »F8 .1 I 6790 

9 .40 FORHAT ( 14X. 1H( . F6. 1. 1H. >F8. 1. 1H) I 6600 

9J9C FURHAT (3(F6. 1.2X11 6610 

ENO *620 
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JJJiginal page is 
of poor quality 


SUIROUTlNR RESULT 


unit is 


fin 6.7*403 


N/M/U. 


1C 


IS 


20 


25 


30 


3S 


CO 


65 


SUBRUUTINE RCSULT ( IT. X.N.F.ANIl, PNALTY, CLOSE, THRESH! tilt 
COMMON NTAAJ«NMAP,NSkG,XH<3,6),YH(3«6I,ARRAYI»76,9I,SP0S1T IS, 100, 36640 
1) , ZC( 3 ) , ZF ( 3 1, NPOS ITS 1 3 ) 6850 
COMMON /PRINT/ SXCENTA(3,3),SYCENTA(3,3 I,SXT1(3,4),SYT1(3,4>,SXT2(6660 


13»4|,SYT2(3,4>,SANSLE<3.3) 

DIMENSION XEN.lt 
00 10 l * l.NTRAJ 
NSEC1 ■ NSEG-1 
00 10 J « 1.NSEG1 

00 .0 K « 1.2 

1 • (l-l)»2»(NSsG-l)*iJ-l)62.K 
IF IK.E0.1I XMI I, J«l) ■ XII. II 
IF (K.E0.2I YHII.J.l) • X(L. 1) 

10 CONTINUE 

WRITE 16.91101 it, f.anh, pnalty, close, thresh 
DO 3U I • l.NTPAJ 

WRITE (6,96201 I, XMI 1, 1 1, YMI 1, 1) 

PII • l./ATANll.i/C. 

00 2C J • 2.NSEG 


6*70 

Aieo 

6190 

6900 

6910 

6920 

6930 

6960 

6950 

6960 

6970 

6960 

6990 

7000 

7010 


20 


WRITE 16.90301 XM(I,J),YM(I, J I. SXCENTR (I»J-l)»5YCSNTR(I»J“l I .$7020 
XTl(l,J),SmiI, J),SKT2(I, J),SYT2(I,Jl.SANGLE(I» J-1I»PII*180. 7030 


30 


7060 
7030 
7060 
7070 
7060 
7090 
7100 
7110 
7120 
7130 

901G FORMAT (1H1.9X.11H ITERATION I .1 2. /. 12X. 12HT0TAI COST? , IP E 16.9, /, 17160 
12X.17HANNOTANCE INIIH . 1PE16.9. /. 12X. 31HPENAITY ON DYNAMIC C0NSTR71S0 
2AINT) >1PE16.9,/,12X.36HPENALTY ON SEPARATING CONSTRAINT! .1PE16. 97160 
3»/»12X,28HPLNAlTY ON THRESHOLD NOISE! .1PC16.9./I 7170 

9020 FORMAT 1 12 X, 16HFLIGHT PATH NO! » 1 1, /. 16X. 10HC0RNER PT., 16X, 16HCENT7180 
lcK OF C1RCLF.6X.15HTANGENTIAL PTS..33X, 13HANGLEIDEGREE ), /,14X, 1H( , 7190 
2F 5. 1, 1H. , . .» 1H 1 1 7200 

9030 FORMAT 1 16X, 61 1H I » FS . 1 . 1H, . F6.1. 1H )• SX I ,F6.1 1 7210 


CONTINUE 

WRITE 16. 9060 1 XM(I. NS EG.ll.YMII.NSEG.il 
WRITE (6.90301 
NPtiSITI • NPOS ITS (I I 

WRITE (6.9060) ( ( SPOS IT ( I, J >K I ,K*1 » 3 1 , J.l.NPOSITI I 
WRITE ( fc.9070) {(SPOSITll.J.KI.K. 1.3). J.l.NPOSITI) 
CONTINUE 

WRITE (9.9069) ( ( ARRAY! I. J I . J*4,5 I . IM.NNAP ) 

RETURN 


9060 FORMAT ( 16X . 1H ( , F8. 1. 1H, »F6. 1. 1H I ) 

9)50 FORMAT 1 / /» 12k, 36HFL I GHT TRAJECTORY X. Y« 
16X.1HY.16X.1HZ./) 

9060 FORMAT I 16X. IPE10. 3, 5X. Ell . 3.3X. E10.3. 5X ) 
9070 FORMAT (3(F8, 1,2X1) 

9)80 FORMAT ( 2 ( 1PE 10. 3, IX ) I 
tNG 


7220 

Z, IN METERS, /.19X.1HX, 17230 
7260 
7250 
7260 
7270 
7280 
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62 


SUBROUTINE CROSOVR 


73/172 r$ 


FTtt 4.7*419 


60/04/24. 10.20.}) 


10 


1) 


20 


2» 


30 


39 


40 


49 


SUBROUTINE CROSOVR (NPOSIT.POSIT.FX) 72)0 

COnnQN KTRA J.NHAP.NSEG, XtH 3, 6), YN( 3,61, ARRAY (976,9 1, SPOS I T( 3, 100 , 37)00 
ll.ZOf }I,ZF'.3),NPCSITS(3I 7310 

CUNHON /CROSS/ X6EGIN.XFINAL,YDIS,CtitITl,CWtIT2 7)20 

DlhlNSIGN CROSS 1 7. 10 1 • POSIT) 3, 100,3 1, NPOSITO) 7330 

Y1NTRP(X,X1,Y1,X2,Y2,X3,Y3«X4,Y4) « < X- X2 1« (X-X)|*(X-X4 1 *V1/ (X1-X27340 
1I/(X,-X)|/(X1”X4)*(X-X1I»<X-X)>*(X-X4I*Y2/(X2>X1I/(X2-X))/(X2-X4I*7)90 
2<X->Xl)*(X-X2>*tX-X4l»V3/(X3-Xl>/tX3-X2l/(X)-X4>*(X>Xl)*(X>X2)*(X-X7)40 
)))*Y4/IX4-Xll/(X4>X2l/(X4-X)t 7)70 

FX * 0. 7)10 

SAHPLE ■ (XBEG1N-XFINALI/11. 7)90 

CROSSU.lt ■ XBf GIN-SAHPLE 7400 

03 10 1 • 2.10 7410 

CROSSU.lt * CROSS tl. I-lt -SAHPLE 7420 

10 CONTINUE . 7430 

00 70 1 * 1.NTR4J 7440 

00 60 J • 1.10 7490 

NP1 . NPCSIT)I)-1 7440 

DO 4u K • l.NPl 7470 

SIGN ■ (POSITU.K.l )-CRQSS( 1, J I !• (POSIT ( I,K*1« 1 l-CROSSU, J 1 17410 



IF (SIGN.GE.g.Ot GO TO 20 

7490 


KK • K 

7900 


GO TO 90 

7910 

20 

IF ( SIGN, GT. 0.0 > GO TO 40 

7920 


IF (POSITd.K. It. NE. CROSSU, Jtt GO TO 30 

7930 


CROSSU* I, Jt > POSIT 1 1.K.2) 

7940 


GO TO 6t 

7930 

30 

CROSS ( 1*1. Jt » POSIT ( I,K*1,2 ) 

7960 


GO TO 6C 

7970 


40 

5j 


63 

70 


99 


CONTINUE 7960 

IF (RK.EQ.lt CRGSS(l*l,Jt • YINTKP( CROSSU. J I . POS IT( 1 . 1.1 1 . P0S7990 
lTU.l»2>»PQSn(I.2.1t.P0SIT(I.2.2t.P0SIT<1.3.1t.P0SlTU*3.2t.76OO 
PUS1T(1«4.1I.P0S1TU.4.2I) 7610 

IF ( KK. NE. 1. AND. KK.NE .NP1 1 CROSSU*I»JI • YINTRP) CR0SSI1. 3 1 . P07620 
SIT ( I.KK-l, 1 ) . PCS1 T ( l.KK-1, 2 1.POSITU.KK.ll.POSIT ( I.KK.2I « POS 1 7630 
1(1,KK*1,1),P031T(I.KK*L,2>,P0$ITU,KK,2,1),P0S1T(I.KK*2,2)) 7640 

IF IKK.E0.NP1) CRUSS(1*I» J) • YINTRP (CROSS U.J )> POS IT II. KK-2, 17690 
t.POSlT(I.KK-2.2t.P0SIT(I.KK~l.ll.PCSIT( I.KK-1. 21. POSIT U.KK. 17660 
I,PUS1T(I,KK,2)»POS1TO,KK*1,1),POSIT(I,KK*1,2)I 7670 

CONTINUE 7660 

CONTINUE 7690 

7700 

7710 

7720 
7730 
7740 
7790 
7760 
7770 
7780 
7790 
7600 
7610 
7620 
7830 
7640 
7630 


NOS TEST THE NEARNESS OF OR THE CROSSOVER BETWEEN TRAJECTORIES 
TYPICAL VALUE CVEIT1 • CWEIT2 • 0.03129 


NACU1 • NTRAJ.l 
00 130 I . 2.NTRAJ 
DO 12u J • 3.NA001 

IP (I.GE.J) GO TO 120 
00 UG K . 1.10 

01S1 • CROSSU, R)-CROSS(J.KI 
IF (ARS(OISI).GE.YOIS) GO TO 80 


FX • FX*evEIU*(Y0IS-DISll**2 


63 


30IRQU7INE IRUSOlR 73/172 TJ 


FIN 1.7*169 


10/01/21. 10.20.13 


69 


69 


70 


6w 


00 


100 

110 

120 

130 


IF (K.ES.IQI GO 70 100 

0132 • CROSS(l»K*lt-CROSS(J.K«l> 

IF (<0131*01321. GE.C. I GO TO 100 
11 * 1 

IF (K.l.ll.CT.lO) GO TO 110 

0133 ■ CRQSS<I»Ktl*Il>-CROSS(J*K«l*ll» 
IF (0152*0133. IE. 0. I GO TO 110 

IF <A8S(0IS2(.IT.A8S<D153H 0132 • 0153 
11 • 11*1 
GO TO VO 
CONTINUE 
GO TO 12b 

FX • FX*CWEIT2*DI53**2 
CONTINUE 
CONTINUE 
RETURN 
END 


7660 

7670 

7660 

7690 

7900 

7910 

7920 

7930 

7910 

7950 

7960 

7970 

7960 

7990 

6000 

6010 

6020 


190006 CH STORAGE USEO 


1.256 SECONDS 
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0P I S? ,AL PAGE IS 
°P POOR QUALriY 


SUAR'JUT (Hi INEwTON 73/172 TS 


PTN A.7*A#5 


• 0 / 0 A/ 2 A. 1 C.I 9 . 3 S 


SUBAOUUNL ONEWTON (HAXIT.STGPCHG, N, KNOW, CELTAA I 


THIS UPTIHIIATION EHPLOTS SELF-SCALING, RESTARTING, 
JUASI-NEWION HETH30. 

REFERENCE) P.C. LUENBERCER INTRO. TO LINEAR AND NONLINEAR 

RRCGRAPPlNG) F. 20A. SEC. 9.3 

H AN 1 T i HAAIPUM NUPBER OF ITERATIONS ALLOHEO 

SIOPC HO STOP IF PERCENTAGE CHANCE IN SUCCESSIVE COSTS IS 

LESS THAN THIS VALUE 

N) OINLNSION CF THE UNKNOWN X 

KNOW) PRESENT OR INITIAL VALUE OF UNKNOWN A 


0 [PENSION ANOHIN. 1). CELTAAINI 

OlNcNSION GNOWOO.l), SNEXTUO.il. POO, l), 000,1), P8I1.1I 
OIPENSICN 6*8(1,11, PPU0»3Ji. 5855(36,30. 5(30,30) 

GIPENSIuN XTEHPUO.l), PTIl.SO. 0SI1.3C), SOU, 301, 0TI1.3OI 
OIHENS ION 5081 30.30 1 
DIHENSION 0(30,11 
IT • 0 
GO TO Ah 
I w FNUV • FSPALL 
00 20 I • l.N 

XNuw (1,1) > ITENPtl.n 
2J CONTINUE 
IT • 1 T * 1 

IF (IT.LT.NAAIT) GO TC AC 
IF (FNALTT.ST.O.) W TO 30 

CALL RtSULI (IT, XNOn.N, FNUW, ANI I, PNALTY, CLOSE, THRESH) 

KeTURN 

30 WRITE (6,9010) PNALTT 
RETURN 

AO CALL CuSTl (0,H,FN0W,AN0W,ANII, PNALTT, CLOSE, THRESH) 

IF (PnALIT.ST.O.) go to so 

CALL POMT (IT, XNQW.N.FNOW, AND, PNALTY, CLOSE, THRESH) 

IF (tUXIT.EO.GI STOP 
GO TO 60 

50 WRITE (6,9w20t FNAETY 


STEP U SeT S • IDENTITY HATRIX AND CALCULATE GRADIENT G 


60 00 TO I • l.N 

CO 70 J • 1,N 
S(1,J) • C. 

IF (UtJ.JI SII.JI • 1. 

70 CONTINUE 

CALL FGkAO (N, KNOW, ANOw, GNUW, OELTAX) 


65 


SUIRJUIINE CNEaIUN 


71/172 


1 $ 


FIN 4. 7*44f 


♦0/04/14. X0.29.lf 


49 


69 


74 


7} 


♦ 9 


89 


90 


99 


140 


1J9 


1 U 


c 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 


. STEP 21 SCI 0 • -SO . 1400 

• . 1610 

1420 

♦410 

♦ 0 cut MPLT (N.N.l.S.GNOW.0.10.30.11 1440 

00 9 k 1 • l.N 1490 

DU. II • -DU. II 1640 

90 CONTINUE 9670 

♦ 440 

4490 

• . 4700 

. >1tP Jl HOST IMPORTANT* . 1719 

. LINE SEARCH ALONG D TO FIND AFA THAT SATISFIES P0*0 . 4720 

. IF COST FUNCTION IS ANALTTIC. SLIGHT MODIFICATION IS . 4710 

. NEEDED (SIMPLIFICATION! (OPTIONAL! . 4740 

• . 4790 

4760 

4770 

K • 0 4740 

ICC K • K*1 4790 

CALL L1NESCH (K, N» PNON. KNOW. 0. AFA2.XTEMP. FSMALL. ANl I. PNAITT. CLOSE, 8800 
ITHHcSHI. RETURNS (10. 1101 4410 

.10 CALL rGFAU (N, FSMALL, XTEMP. GNEXT.DELTAX) 4420 

CD 12C 1 • l.N 4410 

9(1.1! > GKEXT ( 1. 1 )-GNQb( 1.1) 4640 

Mill! • AFA2.0II. 1) 4490 

120 CONTINUE 4860 

CALL IKNSPUS (JO.l.P.PT) 4870 

CALL MPLT (l.N.l.PT.O.PQ. 1.30.1) 4C40 

IF (PU(l.l).GT.C.I GO TO 14U 4890 

IF (K.UT.4) GO TO 10 4900 

PNOR • F SMALL 8910 

DO 131 I • l.N 4920 

XNUNtl.ll • XTEMPU.l) 4930 

13C CONTINUE 4940 

GO TO lWD 8999 

140 DO ISkr t • l.N 8963 

XNUR(,I»1I • XTEMPU.l! 4970 

Hu CONTINUE 4980 

IT ■ 1 T .1 8990 

PRCNT o ABS UF SMALL -F NON ) /FNOW I 9000 

IF IPKCNT.GE.STOPCHG) GO TO 170 9010 

8AI1E (6.4030) PRCNT. STOPCHG 9020 

IF (PRALTT.GT.G.) GO TO 160 9030 

CALL RESULT ( I T , XNUk. N. F SMALL. ANI I .PNALTT. C LOS L. THRESH! 9040 

9090 

RETURN 9060 

160 URITE (6.9C1C) PNALTT 9070 

RETURN 9080 

17 J IF ( II. LI. FAXIT) GO TO 190 9090 

KRITt (O.9U40I 9100 

IF (PNALTT. GT. 0. ) SC TO 180 9110 

CALL ReSULI IK > XKOr.N. F SMALL. ANI 1. PNALTT. CLOSE. THRESH) 9120 

9130 

RETURN 9140 

180 aRi TE (O.9C10) PNALTT 9190 

RETURN 9160 


66 


ORIGINAL PAGE IS 
OF POOR QUALITY 


3UIR3UT1NE CNENIUN 71/172 IS FTN 4.7*499 90/04/24. 10.29.11 


119 141 IF (PNALTY.GT.3.) GO TO 200 *170 

CALL MONIT ( IT. XNOH.N.PSMALL.AMI I.PNALTY.CU3SE. THRESH! 9160 

GO TO 210 *190 

2(0 WRIT i (O.9C20I PNALTY 9200 

210 CONTINUE 9210 

120 C 9220 

C i ...................... 9210 

C . . 9240 

C . STEF 41 IF FIT* IS INTEGER MULTIPLE OF N GO TO STEF ll • 9290 

C . IF NUT. UPDATE S , 9260 

129 C . IS »IT* INTEGER MULTIPLE OF K4 . 9270 

C . . 9290 

C 9290 

C 9100 

IF l(fLC)ATlIT)/FLO»T(NII .NE. FLOAT ( IT/N) ) GO TO 220 9110 

llw FNOW > FSMALL 9120 

GJ TO 6o >330 

C 9340 

C 9390 

C . . 9360 

139 C . UPDATE MATRIX S) GO TO STEP 2 . 9170 

C . . 9360 

C 9390 

C 9400 

220 CALL MPLY IN.N, 1, S.W, SO, 30,30, 1 ) 9410 

140 CALL 1RNSP0S (30. 1.0, ST) 9420 

CALL MPLY (N.l.N.SO.GT.SQQ, 30.1.30) 9430 

CALL MPLY (N,M,N>S00>S,SaaS»30.3fi.33) 9440 

CALL MPLY (1,N,N,OT,5,OS,1,39,3G) 9490 

CALL MPLY (l.N.l.OS. 0.0S0(l,3C.l) 9460 

149 CALL MPLY (N.1.N.P.PT.PP.30.1.H-) 9470 

00 23w 1 • l.N 9460 

DO 23U J • l.N 9490 

sd.j) • (S<i,J>-saas<i,j>/oso<i,n i* (pot i,i > /osQ< i.i >).pp<i, j 1/9900 
A PGU..I 9910 

190 230 CONTINUE 9920 

FNOW • FSMALL 9930 

00 240 I • l.N 9940 

GNUM 1,11 • GNEXT (1,1) 9990 

240 CONTINUE 9960 

199 C-0 TU SC 9970 

C 9980 

9«il</ FORMAT ( 9X.43HIN RESULT! DYNAMIC CONSTRAINTS VIOLAT ION, / , 10X.9HPEN9990 

1..LTY •» 1PE16.9) 9600 

9020 FORMAT (9X.39HIN MONITl DYNAMIC CONSTRAINTS V ICLAT10N, /, 10X, 9HPENA9610 

160 an TY ..IP616.9) 9620 

9430 FORMAT (2X.37HPF RCENTAGE CHANCE IN SUCCESSIVE COSTS, 1PE10.3.29H LE9630 

ISS THAN STOP CRITERION. 1PE10.3I 9640 

9340 FD.MAT (2X, 29HMAX 1MUM ITERATION SET REACHED ) ' 9690 

END 9660 

41U0W3 CM STORAGE USED 1.246 SECONDS 
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j 
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'*T ( 


SUBROUTINE N?L» 73/172 TS 


FTN A.7M89 BO/OW2A. 10.20.SS 


J 


» 


1 


S 


10 


IS 


SUBROUTINE mr (L » M« N> A> 8»C« LDEC >HDEC> NOEC ) 9670 

C 9660 

C 9690 

C . . 9700 

C . CALCULATE MATRIX MULTIPLICATION C • A« . 9710 

C . . 9720 

C 9730 

C 97A0 

DIMENSION ALLCEC.HDEC). BLMDEC.NOEO# C(LOEC.NDEC) 9790 

DO 10 I • 1.. 9760 

DO 1U J • l.N 9770 

C(I.J) • 0. 9760 

DO IU K • 1,M 9790 

LII.J) • C(Ii>J) + A<I»K)»BtK.J) 9600 

10 CONTINUE 9610 

RETURN 9620 

END 9630 


AluJCB CM SICRAuL JSED 


. IAS SECONDS 
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l 


SUBROUTINE TRNSPUS 


73/172 


TS 


PIN 4. 7*4(5 


•0/04/24. 10.20.15 


SUBROUTINE TRNSPOS IHfNfAfBI 9*40 

C 4(50 

C 

C . • 4970 

5 C . TRANSPOSE CP PURI* 4 IS RETURNED IN NATRI* S • 9B(0 

C . • 9B90 

C 9900 

C 9910 

DINENSION A( N. N ). BtNf HI 9920 

10 DU 10 I • 1. N 9910 

DO 10 J • It N 9940 

B(Jfl) • AII.JI 9950 

10 CONUNUE 9940 

RETURN ”7 0 

15 END 4980 

41000B CN STORAGE USED .090 SECONDS 
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ORIGINAL PAGE 19 
OP POOR QUALITY 


* i 


* 


f 


SUBROUTINE f 5(1*0 73/172 TS 


FTN *.7«*69 10/0*72*. 10.20.33 


-f 


,y- 


ts 


SUBROUTINE FGRAD (N. F» X.G. BELT** I 9990 

C 10000 

C 10010 

C . . 10020 

9 C . CALCULATE GRADIENT OF COST F WITH RESPECT TO UNKNOWN X . 10030 

C . • 100*0 

C 10090 

C 10060 

DihENSIUN X(N# 1 )# Gt30.ll. DELTAXi 301 10070 

io on io i ■ i»n looio 

XU. II ■ X ( 1. 1 (.DELTAXt 1) 10090 

CALL C0ST1 (l.N.FF.X. ANI I.PNALTY.CLOSE. THRESH) 10100 

6(1.1) • (FF-FI/BELTAXtl) 10110 

X(I.l) • XII. l)-DELTAX(l) 10120 

19 10 CONTINUE 10130 

RETURN 101*0 

END 10190 

*16066 CN STORAGE USED .137 SECONDS 


1 

I 

I 

1 

1 

I 

I 

I 

1 

1 
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J 

I 

J 

J 


SUBROUTINE GAFA 73/172 TS 


FIN 4.7*469 10/04/24. 10.20.19 


i 

I 

I 


9 


10 


19 


XUBROUT IKE GAFA I N, F, OAFA,G, X,D I 
C 

c , 

c . 

C , CALCULATE GRADIENT OF F WITH RESPECT ?G AFA 
C . 


C 

01 HENS ION AIN, II, 0130,11# ATENPI30,1I 
OU H I ■ 1, N 

ATEHF 11,11 • X(I,1)«DAFA»D(I,1) 

10 CONTINUE 

CALL C0ST1 tl,N,FF,XTEHP«l,ANIl,PNALTV,CLOSE,FXI 

G ■ (FF-FI/DAFA 

RETURN 

END 


101*0 
10170 
. 10110 
. 10190 
. 10200 
• 10210 
. 10220 
10230 
10240 
10290 
10260 
10270 
10260 
10290 
10300 
10310 


i 


410006 CN STORAGE USED .104 SECONDS 


I 


I 

l 
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J 


SUBROUTINE EMHOR 73/172 TS 


FTN 4.7«4*» *0/04/2*. 10.20.1* 


i 


SUBROUTINE ERROR 
WRITE (6,10) R 

stop 


(K) 10320 

10310 
10340 
103*0 

10 FOR HA 7 (1X,6HAFTER ,U,33H TIMES THROUGH LINE SEARCH, STUl,37H CA1Q160 
1NNQT FIND AFA WHICH SATISFIES PK>0) 10170 

END 103*0 


41000* CM STORAGE USEO 


.036 SECONOS 
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SUBROUTINE LINESCH 73/172 7S 


FTN «.T«U9 10/04/24. 10.29.19 


5 


10 


IS 


20 


2S 


30 


35 


40 


45 


9b 


55 


SUBROUTINE LINESCH ( K, N, FNOW, XNOW, 0, AFA 2, XTENP. FSNALI.ANII.PNAITY. 10390 


1CLQSX. THRESH). RETURNS (N.NNI 10400 

01NENS1QN U10.lt. DI30.lt. XTENM30.lt. XNCWIN.lt 10410 

OINENSICM XK30.lt. X2(30.11> XM30.lt. X4dO.ll. AFA(4I» FXdt 10420 

C 10430 

C 10440 

C . . 10490 

C . CUBIC FIT BY INITtXLLY LETT INC 4FA • 0 AND DAFA • .01 . 10460 

C . . 10470 

C 10410 

C 10490 

IF (K.LT.-lCbi GO TC 60 10500 

ONAX • 0(1.11 10510 

00 10 I • l.N 10520 

IF (DHAX.LT.A6S(0(I.l)t) ONAX • ABS(0(I.ltt 10330 

10 CONTINUE 10540 

IF (UNAX. £b. 0.1 CALL CHECK ( l.N. FNOW. XNOW, 0. AFAO.FAFAO.GAFAO.AFAl. 1055?* 
lFAFAl.GAFAl.Ul.UUl.U2i 10560 

OAFA • 10. /ONAX 10570 

OAFA • b.Cl 10360 

AFAb • 0. 10590 

FAFA9 • FNOW 10600 

CAUL GAFA IN.FAFAO.OAFA.GAFAQ.XNOW.O) 10610 

AFA1 • ICu./ONAX 10620 

00 2b I • l.N 10630 

XVENFll.lt • XNOW ( I. 1 llAFAl*D( I. It 10640 

2w CONTINUE 10650 

CALL C0ST1 (0.N.FAFA1.XTENF.ANII.PNALTY.CLCSE.THRESH) 10660 

3b CALL GAFA (N. F AFA1.0AF A.GAFA1. XTENP. 01 10670 

IF (ABS((AFA0-AFAi|/AFAlt.LT..001l CALL CHECK (2.N.FN0W. XNOW.O.AFA106BO 

Ib.FAFAb.GAFAC.AFAl.FAFAl.GlFAl.Ul.UUl.U2t 1069b 

U1 • GAFA0»GAFA1-3.«(FAFAC-FAFA1I/(AFA0-AFA1) 1070C 

UUi • U1»*2-GAFA0*GAFA1 10710 

IF (UU1.LE.6.) U2 • 0. 10720 

IF (UUa.GT . 0. ) U2 • SORTlUUlt 10730 

IF (GAFA1-GAFAO*2.»U2.EJ.O. ) CALL CHECK d. N. FNOW. XNOW. 0. AFAO. FAFA10740 
iJ.GAFAC.AFAl.FAFAl.GAFAl.Ul.UUl.U2) 10750 

AFA2 • AFa1-(AFA1-AFA0)*(GAFA1«U2-U1)/(GAFA1-GAFA0»2.*U2) 10760 

DU 4w I • i.N 10770 

XTENP ( I. 1 ) • XNOW ( I. l)*AFA2*D( I . 1 > 10760 

4b CONTINUE 10790 

CALL C0ST1 (C.N.FSNALL, XTENP, ANII.PNAITV, CLOSE, THRESH! 10800 

IF ( FSNALL .GE.FNOW ) GO TQ SO 10S1C 

CALL CUSTi ( 1, N. FSNALL. XTENP, AN I l.PNALTY. CLOSE. THRESH) 10B20 

RETURN NN 10830 

5b AFAu • AFA1 10640 

FAFAb • FAFA1 10890 

GAFAb • GAFA 1 10869 

AFA 1 • AFA2 10870 

FAFA1 • FSNALL 10880 

IF (APA2.E0.0. t DAFA • b.Ol 10890 

IF (AFA2.NE.0.) DAFA • 0.G1*AFA2 10900 

GO TO 5 J 10910 

C 10920 

C 10930 

C . • 10940 

C . CQGGIN AL- > (HN TOUACRATIC FIT) . 10950 


73 


SUttHOUT 1H£ UNESCH 73/112 


TS 


FTN 4.7*415 


(0/04/24. 10.20.15 


C . . 10450 

i 10970 

01 C 109(0 

tv DM41 • 011.11 10990 

DO 7C I • l.N 11000 

£1(1,1) • XNOWII.l) 11010 

IF (0HAK.L7. 011.11) OHAX • 01!. II 11020 

65 7w CONTINUE 11030 

04F» • SOw./QMAX 11040 

FI • FNOK 11090 

SION • 1. 11060 

SUP • 1. 11070 

Ttl 00 bG 1 • l.N 110(0 

X2I1.1) « XH5U! «.1)*D4F 6*0(1. II 11090 

bO CONTIlH.-.. moo 

C*U COS?i f! i l 1**72. 4N11.PN4LTY, CLOSE. tHRtSHI 11110 

9t IF IF2. (*.>».', ,* U 110 11120 

75 STEP . 2, 11110 

DO 100 1 • 11140 

£111,11 - £2*1.1) 11150 

X2I1.1I • £211. 11*2. *D4FA*D<I»1I*S1GN 11160 

100 CONTINUE 11170 

(tl FI • F2 11160 

CALL C0ST1 (1.N.F2.X2.AN11.PNALTT. CLOSE. THRESH) 11190 

GO TO 40 11200 

110 IF (STbP.NE.l.l GO TO 130 11210 

STEP • 2. 11220 

(9 SIGN • -1. 11230 

DO 120 1 > l.N 11240 

X2I1.1) • XNCNI1.1)*DAFA*D( J.1I«S1GN 11290 

1 20 CONTINUE 11260 

CALL CUST1 II. N.F2.X2.AN11.PNALTY, CLOSE, THRESH) 11270 

90 GO TO Vu 112(0 

130 DO 140 1 • l.N 11290 

K 3 ( 1,1) • 1X1(1, l)*X2(1, 1)1/2. 11300 

140 CONTINUE 11310 

CALL COSTl 11, N.F3.X3.ANI1.PNALTY, CLOSE, THRESH) 11320 

99 AFAllI • 1X1(1, 1I-XNDWI1, 111/0(1.11 11330 

AFAI2) • 1X2(1, ll-XNOWU. 111/0(1. II 11340 

AFAI3I • (X3U, ll-Xf.QkU.il 1/0(1,11 11350 

ISO AFAI-,1 • 0.S*(((AFA(2I**2-AFA(3I**2I*F1*(AFA(3I **2-AFA ( 11**2 l*F2*( 11360 
lAFA(l)**2-AFAt2|**2l*F3)/((AFA(2l— AFA(3I)*F1*(AFA(3I-AFA(1))*F2*(A11370 
110 2FAU)-AFA(2))*F3I) 11380 

DU 160 1 • 1,N 11390 

X4(l,l) » XN0kU,l>*AFA|4)»D(I,l) 11400 

16C CUNUNUE 11410 

CALL COSTl II, N,F4»X4»ANI1»PNALTY» CLOSE, THRESH) 11420 

115 FXI1) • FI 11430 

F X 1 2 ) • F2 11440 

FXI3I • F3 11490 

FM1N • FX(l) 11460 

HIN • 1 11470 

lit, MAX • FX 1 1 1 11480 

MAX • 1 11490 

OJ id*. I • 2,3 11500 

IF (FMIN.LT.FXU) I GO TO 17w 11910 

FM1N • FX 1 1 1 11520 


SUIROOI INt 

LINESCH 73/172 TS 

FTN 6.7*663 

•0/06/26 

113 


KIN • 1 


11330 


170 

IF If fUX.6T.FXU I 1 60 TO 166 


11360 



FIUX • mil 


11330 



fUX • I 


11360 


160 

CONTINUE 


11370 

120 


DO 190 I • l.N 


11360 



XCCINF • XNQWII»1>«AFA<N1N)*D<1»1I 


11390 



IF IABS ( XC0HP-X6I I.1II.IE.1C.I 60 TO 200 


11600 



XTENPIl.il • XCOHP 


11610 


19 C 

CONTINUE 


11620 

125 


CALL C0S11 IO.N.FSNAIL.XTENP. ANI I. PNALTT.CLOSE. THRESH 1 

11630 



RETURN NN 


11660 


200 

AFAIHAXI ■ AFAI6I 


11630 


GO TO Hi H660 

fcND 11670 


610 106 CH STORAGE USED 


2.336 SECONDS 


I 


.4 


SUBROUTINE CHECK 73/172 15 


FTN 6,7*685 SO/C6/26. 10.20.3? 


5 


10 


IS 


SUBROUTINE CHECK ( I CHECK; N. FNOW. XNOW.D, AFAO.FAFAO.GAFAO, AFA1.FAFA111680 
1,GAFA1,U1,UU1,U2) 11690 

DIMENSION XNOW(N.l). 0(30.1) 11700 

WRIIc (6*101 FNOV.N,(I»XNOV(l,l).D(I»l)»I*l.N) 11710 

IF ( ICHcCK.EQ. 1 ) STOP 11720 

WRITE (6.201 AFAO.FAFAU.GAFAO.AFAl.FAFAl.GAFAl 11730 

IF ( ICHECK.E0.2 ) STOP 11760 

WAITE (6.301 U1.UU1.U2 117S0 

STOP 11760 

11770 

10 FORMAT (20X.66H1HIS IS SUBROUTINE CHECK WHICH GIVES ALL THE.3SH INU780 
IFOR MA HON IN SUBROUTINE L INESCH. . / /. 33K.7HFN0W • . IPE16. 9. //.3TX, 611790 
2HXNUW,17X»1HD,//»*(29X.12,3X.1PE16.9,6X,1PE16.9 ) ) 11800 

20 FORMAT ( / /. 20X. 7HAFA0 • .1PE16.9.9X.6HFAFA0 • > 1PE16.9. 5X. 8HGAFA0 11810 
I* .1PE16.9./.20X.7HAFA1 • . 1PE16. 9, 5X, 8HFAFA1 ■ * 1PE 16.9, SX. BHGAFAll 820 
21 • .1PE16.9) 11830 

30 FORMA I ( //.20X.5HU1 • . IE16.9.5X.6HUU1 • . 1PE16. 9, 9X, SHU2 • .1PE1611B60 
1.9) 118S0 

t N0 11660 


610008 CM STORAGE USED 


.163 SECONDS 


* 
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